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Abstract 

The electromagnetic (EM) form fators of the pseudoscalar and vector mesons 
are calculated in the light-front framework. We extract these form factors 
from the relevant matrix elements directly, instead of choosing the Breit frame. 
The results show that the charge radius of the meson are related to both the 
first and second longitudinal momentum square derivative of the momentum 
distribution function. In addition, the static properties of the EM form factors 
and the heavy quark symmetry of the charge radii are checked analytically 
when we take the heavy quark limit. 
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The understanding of the electromagnetic (EM) properties of hadrons is an important 
topic, and the EM form factors which are calculated within the non-petrturbative methods 
are the useful tool for this purpose. There have been numerous experimental and 
theoretical studies [0-0 of the EM form factors of the light pseudoscalar meson (vr ana K). 
However, duing to the experiments are more difficult, the EM form factors of the light vector 



meson (p and K*) have fewer investigations than their pseudoscalar counterparts ||TT| , p!2 
even though they could provide much information about the bound-state dynamics. As for 
the EM form factors of the heavy mesons (which containing one heavy quark), there are 
much fewer study than the light ones. In the heavy hadron investigation, however, the heavy 
quark symmetry (HQS) [0 is a fundamental and model-independent property. In this work, 
we will study the EM form factors of the heavy mesons with the light-front framework, and 
will check whether HQS is satisified or not among these EM properties. 

As is well known [|14| , the light-front quark model (LFQM) allows an exact separation in 



momentum space between the center-of-mass motion and intrisic wave functions, a consistent 
treatment of quark spins and the center-of-mass motion can also be carried out. Thus it 
has been applied in the past to calculate various form factors [|T5|. Let us begin with the 
light-front formalism. A meson bound state consisting of a heavy quark qi and an antiquark 
q2 with total momentum P and spin 5* can be written as 



|M(P, S, S,)) = I {d^pi}{d'p2} 2{2^f5\P - - P2] 



X ^^'^^(pi,P2,Ai,A2) |gi(pi,Ai)g2(P2,A2)), (1) 

Ai,A2 



where pi and p2 are the on-mass-shell light-front momenta, 

,2 I ^2 



P = iP'^,P±) , P± = (p\p^) , P = ^ t^^ ' (2) 

p+ 



and 



3 dp+(fp^ 

\q{pi, \iMP2, \2)) = &Upi)4.(P2)|0), (3) 
{bAp'),b{{p)} = {d,,{p'),d{{p)} = 2{2'nf 6%p' - p) 6,,,. 

In terms of the light-front relative momentum variables (x, k±) defined by 

p+ = (l-x)P+, p+ = xP^, 

Pi± = (1 - x)P± + k±, p2i_ = xP_L - k±, (4) 

the momentum-space wave-function vE''^'^^ can be expressed as 

vl>^^^(pi,P2, Ai, A2) = RfXix, k^) 0(x, A;^), (5) 

where (f){x, k±) describes the momentum distribution of the constituents in the bound state, 
and -Rff^j constructs a state of definite spin {S, S^) out of hght-front helicity (Ai, A2) eigen- 
states. Explicitly, 
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Rx%{x,k±) = J2 (Al|7^j^^(l -a;,/cx,ml)|sl)(A2|7^j^^(x, -k±,m2)\s2){-Si; -S2\S,S^), (6) 



Sl,S2 



where \si) are the usual Pauh spinors, and TZm is the Melosh transformation operator |]17 

rrii + xMq + ia ■ k± X n 



TlMix,k±,mi) 

^{rrii + xMo)2 + kl 

with n = (0, 0, 1), a unit vector in the 2;-direction, and 



1-x) 



X 



In practice it is more convenient to use the covariant form for -Rff^^ [EH- 



P1P2 



V2 Mo 



u{pi, \i)Tv{p2, A2), 



where 



Mo = ^jMi - (mi - m2)2, 

r = 75 (pseudoscalar, 5* = 0). 

^ • {Pi -P2) 



V = - ^(5,) + 



Mo + mi + m2 



with 



P 



-£^(±1) ■ Px, 0, £l(±l) 



(vector, 5* = 1), 



1 



-M2 + P 



(7) 



(8) 



(9) 



(10) 



Mo V ^+ 
We normahze the meson state as 

(M(P', S\ S',)\MiP, S,)) = 2{2TxfP+5\P' - P)6s's5sis. , 

so that the normahzation condition of the momentum distribution function can be obtained 

dx (Pk 



(11) 



/ 



2(27r)= 



(x, k^)\' = l. 



(12) 



In principle, (f){x, k±) can be obtained by solving the light-front QCD bound state equation 
p!8| . However, such first-principles solutions have not been done yet, and the phenomeno- 



logical wave functions must be considered. In this paper, we don't designate any one to use 
but only ask the momentum distribution function has the scaling behavior [19 

(pQqix^kx)^ y/m^^{mQX,k±), (13) 
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in the infinite quark mass limit mq oo. The factor y/rrTq comes from the particular 
normalization we have assumed for the physical state in Eqs. ([Tl| ) and (0). The reason 
why the light-front heavy-meson wave function should have such an asymptotic form is as 
follows. Since x is the longitudinal momentum fraction carried by the light antiquark, the 
meson wave function should be sharply peaked near x ~ Aqqb/^q- It is thus clear that 
only terms of the form"mQx" survive in the wave function as itlq —>■ oo; that is, rrigx = X 
is independent of rriQ in the rriQ oo limit. 

Now we consider the pseudoscalar meson first. The EM form factor of a meson P is 
determined by the scattering of one virtual photon and one meson. Since the momentum of 
the virtual photon q is space-like, it is always possible to orient the axes in such a manner 
that = 0. Thus the EM form factor of a pseudoscalar meson P is determined by the 
matrix element 



{PiP')\J+\PiP)) = e Fp{Q^){P + P 



(14) 



where = qegC'-y^q is the vector current, is the charge of quark q in e unit, and = 
(P' — Py < 0. This process is illustrated in Fig.l. With LFQM, Fp can be extracted by 
Eq. (m 



FpiQ') 



dx(fk± 4>p{x,kj_) (f)p'{x,k'j_ 



'J2 



+ e, 
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dxd'^k± (j)p{x,kj_) (f)p'{x,k'lj 



+ k^ ■ k\ 



A^ + k^- k'[ 



(15) 



where k'j_ = k± + xQ±, k'^ = k± — {1 — x)Q±, A = xnii + (1 — x)m2, and J A^ + k 



Mq^x{1 — x). From Eqs. (P), (|^), and @, it is understandable that the term \J A^ + k^ 
comes from the Melosh transformation. After fixing the parameters which appear in the 
wave fnction, Eq. (^) fits the experimental data well |20]. But this is not the whole story. 
We consider the term (pp = (j)p{x, k±)/ J A^ + k\ and take the Tayor expansion around k"]^ 



cj)p,{k'l) = (j)piki) + 



dki 



{k'l - ki) + 



d\ 



2{dk 



2 \2 



{k 
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,2 \2 



+ 



Then, by using the idenity 

j d^k^ (fcx ■ Qx)(A:x ■Q±) = lJ d'k± kl Q] 
we can rewrite fllSD to 



(16) 



(17) 



Fp{Q' 



dxd'^kj 
2(27r)3 



■0p(x, k_i){eq^x'^ + 6^2(1 - x) 



2kl 




+ Qpk 



(18) 



where 
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e 



M 



'M 
2 \2 



(19) 



The mean square radius {r^)p = 6{dFp{Q^) / dQ^)\Q2=Q can be obtained easily that (r^)p 
{r'^)qi + (^^)?2> where 



92 



'92 



-6 
-6 



dxd'^kx 
2(27r)3 

dxd'^k± 
2(27r)3 



(^2 + 2k] 



_d_ 



dki 



+ (^2 + 



dB 



(20) 
(21) 



It is worthwhile to mention that, first, the static property -Fp(O) = ep is quite easily checked 
in ([l8|). Secondly, from Eqs. (pO]) and (^T]) , we find that the mean square radius is related 
to the first and second longitudinal momentum square derivatives of which contain the 
Melosh transformation effect. When we consider the heavy meson qi Q, we get 



X 



X 



A^ A = X + 1712, 



(22) 



and Eq.(| 



in the heavy quark limit (mg 
f -6 r dXd^k 



0, 



12 



-92 



2(27r)3 

dXd'^k_ 
2(27r)3 



oo), then the mean square radius 

2 



(^' + 2h!\ 



_d_ 



+ {A^ + A;1)A;^ 



(^2 + 2A;^ 



+ {A^ + fci)fc^ 



dB, 



(23) 
(24) 



where $ = ^/y^^ + /c^.- Eq. (^) means that, in the heavy quark limit, the mean square 
radius {r'^)p is blind to the flavor of Q. This is the so-called flavor symmetry. 

Next, the vector meson case is considered. As in the pseudoscalar case, we set = 0, 
then the EM form factor of a vector meson V is parametrized by the matrix element 



e e 



{V'{P',e')\J+\V{P,e)) = GriQ')e ■ e'{P + Py 

-G2{Q'){P-e' e+ + P' 

+ G'3(g')7Tl^P-e' P'-eiP + P'y, 



2Ml 



(25) 



where My is the V-meson mass and e(e') is the polorization of ViV'). In the literature 
2T|^4|, these invariant form factors are usually obtained by choosing the Breit frame. For 



the vector meson there are only three independent form factors, whereas four matrix elements 
of are independent since is (i) Hermitian, (ii) invariant under rotations about 
ra, (iii) time reveral, and (iv) reflection on the plane perpendicular to n. An additional 
constraint, which called the angular condition, is related to the rotational invariance of the 
charge density [0. This invariance is related to the unitary transformations based upon a 
subset of Poincare generators. The authors of Ref. p3[ had written the angular condition 
in the following form 
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A{Q') = (1 + 2r/) Ji,i + Ji,_i - JSr^Ji.o - Jo,o = 



(26) 



where t] = Q^/AMy and J^/^e = (^'(-PS ^Ol-^^l^l-^) ^))- generally, this condition must 
be satisfied when the matrix elements of the many-body current are existence, and are 
breakdown for the one-body one alone. Thus, when only the one-body current operator 
g7^g is considered, Eq.(^) usually doesn't equal to zero and the calculation of these form 
factors depends upon the description used |T^. Following p^ , |25[| , Jo,o is expected to be 
the dominant matrix element in the pQCD regime. Thus, as the arguement in [^, it is 
reasonable that we calculate the Gj form factor without using the matrix element Jq.o in the 
constituent quark model. 

In this paper, we will extract the form factors directly, instead of choosing the Breit 
frame. The dependence of these form factors must be calculated order by order. In spite 
of the fact that we cannot obtained the exact form, the advantages, however, are fruitful 
and will be mentioned later. We choose the matrix elements Ji_i, Ji,_i, Ji,o and Eq. (|lj 
to extract the form factors Gi, G2, and G3 in the suitable order 



G,iQ' 



+ e. 



12, 



G2{Q' 



Ml 



V 



dx(Pk_\ 
2(27r)3 



V 1 ^qi 



x'eviA^ + 2fci) + x'QvkliA^ + kl) 



M2 



+ 



+ 0{Q' 



dxd'^kj 
2(2^ 



xa + x^\h + 



2(1 
d 



x)M^W^ 
xk"^ ^ 



+ e 
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W^ 



V . 



G,{Q^ 
AMI 



+ O(Q^), 
dxd'^kx 
2(27r)3 

+ 0{Q'), 



+ xk'^:Qv\b + 



[e^^x^ + 652(1 - xf)(l)y< 



2Wv \ 

+ 



C + 



d 



A 



l-x)MoWv, 
X ^ {1 — x) 

Ak^^ 



Wv 2x(l - x)MoWj 



(27) 



(28) 



(29) 



where W^ 



V 

1 

Mo 
-4 
Mo 



Mq + nil + ^2 and 



1 — X 

X 



1 — X 
nil — 



ml + kl 
ml + kl) - 



X 



X 



m'2 + ki) + M^{l-2x) 
mi + m2 



X 



-2 /_2 I 2 



[ml + ml)) 



{M^ - (mi + mi)) 



d = 4- 



Mo 

X 



Mr 







Wv 
{mi + m2] 



Mo ml + ml 
^ MqWv 



X 



{ml + kl) 



1 



X 



1 — a; ■ * ' X 

r2 I 



ml + k^\ 



+ (Mo - mi - m2){M^ - {m{ + m^)) 



(30) 
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These Gi{Q'^) form factors are related to the charge, magnetic and quadrupole form factors 
as 
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Gem 

Gm{Q^) 
Gq{Q') 



Gi -\- K-Gn, 

—G2, 
G\ + G2 



1 + 



4M2 



\G.. 



(31) 



These form factors are proportional to the static quantities of charge e, magnetic moment 
/i, and quadrupole moment Q when = 0, 



^^(0) = ev, 

Gm{0) = /i, (the unit is e/(2My)), 
Gq(0) = Q M^. 

The mean square radius of the vector meson (r^)y = 6{dGc{Q'^) / dQ'^)\Q2^Q is 



(32) 



dxd'^k I 



6 



+ CqJ - 6 



2(27r)= 



dxd'^kj 
2(27r)3 



(^2 + 2A;2 



+ 



(xW^v)^ ' 2x(l-x)MoW^ 



>v 



:i-x) 



{A^ + 2k] 



_d_ 



kl 



+ 



{{l-x)Wvf 2x{l -x)M^W^ 



H — 775" (^1 + G2 + G3 



(33) 



Q2=o 



Similar to the pseudoscalar case, in addition to Eqs. (|T3D and (^21) , there are some behaviors 
existing 



MQ-^mQ^Oirnqx), Wy ^ 2mQ + 0{mQx) 



(34) 



if the heavy quark limit is taken. We also obtain the Eqs. ( ^3]) and (p4|) for vector meson 
case. Thus the mean square radii of pseudoscalar and vector meson are equal {r'^)p = {f'^)v- 
This means that they satisfy the spin symmetry. Combine this relation and Eq . (p3D , we 
find that the light degrees of freedom are blind to the flavor and spin orientation of the 
heavy quark. This is the so-called heavy quark symmetry (HQS). Up to now, we have not 
designated any wave function yet, this also satisfies the well-known property that HQS is 
model- independent. Reviewing the processes, we can realize that, in this approach, the 
static properties of the EM form factors and the heavy quark symmetry of the mean square 
radii can be checked much more easily than in the Breit frame. This is the major reason 
why we calculate the dependence of those form factors order by order. 
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In conclusion, we have calculated the EM form factors of the pseudoscalar and vector 
mesons. Specially, the EM form factors of vector meson are extracted from the relevant 
matrix elements directly, instead of choosing the Breit frame. We found that the charge 
radius is related to both the first and second longitudinal momentum square derivative of 
the momentum distribution function. We also found the static properties of the EM form 
factors and the heavy quark symmetry of the mean square radii are checked analytically by 
evaluating the dependence of those form factors order by order. Therefore, in the heavy 
quark limit, the charge radii of pseudoscalar and vector mesons have both spin and flavor 
symmetries, and these properties are model-independent. 
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FIGURE CAPTIONS 
Fig. 1 Diagrams contributing to the EM form factors of meson. 



10 



